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Abstract 



In this work radiative corrections in the total hadronic decay rate 
of the r lepton and some moments of its differential distributions are 
studied employing perturbative QCD and the operator product ex- 
pansion. We calculate quadratic quark mass corrections in the strange 
mass to the decay rate ratio R-r to the order C(a^m^) and find that 
they contribute appreciably to the Cabibbo suppressed decay modes 
of the r-lepton. Using the results of a recent experimental analysis, 
we obtain m,(l GeV) = 200 ± 40e^p ± 30th MeV. 



^On leave from Institute for Nuclear Research of the Russian Academy of Sciences, Moscow, 
117312, Russia. 



1 Introduction 



With an ever increasing number of r leptons observed by the four LEP experi- 
ments and by CLEO the Cabibbo suppressed r decays have become one of the 
important topics of recent experimental analysis The study of exclusive chan- 
nels allows to determine hadronic resonance parameters, to test the predictions 
in the chiral limit and, quite generally, to explore the hadronic current in the 
low energy region. Multi-differential angular distributions of mesons can be used 
to measure the polarization of the r lepton and to furthermore disentangle the 
various spin parity contributions of hadronic states with = , 0' , 1" or 1+ 
induced by the (non-) conserved parts of vector and axial vector current respec- 
tively 0. Complementary studies are based on the inclusive decay rate derived 
from the semileptonic branching ratio or the r life time (see e.g. P, ^, |], H, ^). 
The determination of the strong coupling constant a.^ has been a focal point in 
these investigations. The small reduction of the Cabibbo suppressed rate (relative 
to the massless prediction) has been used recently even to extract a value for 
the strange quark mass. The analysis was based on the total rate and a theoret- 
ical calculation ^ including quadratic mass terms with coefficients of order a^. 
Recently, the authors of [§] have revised their analysis and corrected an error in 
the published numerical value of the term. As a result the contribution from 



that term has increased dramatically by a factor of about four pO| , |TT|. This new 
development calls for a fresh look at the possibility for an determination from 
r decays. 

In this work the theoretical analysis will be extended and improved in various 
ways: 

• QCD corrections to the mass terms of order will be calculated for some 
of the moments of the spectral functions which if included lead to a sizeable 
shift of the predictions. 

• It is proposed to use the techniques of [0 to separate the states of dif- 
ferent spin parity, allowing thus for quark mass determinations from four 
independent spectral functions. Eventually even different moments of their 
respective spectral functions might be considered, leading to additional tests 
of the method. 

• The method of resummation of effects from the running of both the 
coupling constant and the strange quark mass along the contour of integra- 
tion in the complex plane through the renormalization group improvement 
is used for r lepton observables to provide better convergence of the per- 
turbative series. 
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The configurations with 0"^ and 0~ and non-vanishing invariant masses are 
strictly forbidden in the limit of massless quarks. The integrated rate, i.e. the low- 
est moment of the spin separated distribution exhibits a remarkable dependence 
on the perturbative mass of the strange quark and a nonperturbative parameter 
at the same time. Higher moments are free from this nonperturbative constant, 
allowing thus, at least in principle, a new determination of the strange quark 
mass. 

The paper is organised as follows. In the next section some general relations 
are given and our notation is fixed. The observables to be analysed in r lepton 
physics are introduced and the stage is set for their perturbative analysis. In Sect. 
3 some new features are described which appear due to the mass terms in the 
Cabibbo suppressed channel. Explicit expressions for the coefficients of quadratic 
mass terms are given. In the same section both the finite order and resumed 
observables are presented. Sect. 4 concentrates on nonperturbative corrections to 
spin zero contribution. In sect. 5 numerical results are given. The last section 
6 contains our conclusion about the possibility and accuracy of a strange quark 
mass determination from an inclusive analysis of Cabibbo suppressed modes. 



2 Generalities 

In order to set the framework of the subsequent discussion, which includes quark 
mass effects and the separation of the spin one and spin zero contributions, we 
repeat below a number of generalities about the theoretical analysis of r lepton 
observables both in massless and massive case and also introduce our notations for 
necessary quantities which can be found in earlier literature (see, e.g. Refs. |1^). 



2.1 Correlators 

Physical r lepton observables are related to correlators of vector and axial vector 
currents of light quarks that are defined as follows 



with = Y^f=u,d,s^} = qi'^nilblQi- Here Qi and Qj are two (generically 

different) quarks with masses and irij respectively. In the present paper we 
work within QCD with effective three light quarks and do not consider corrections 
due to heavy quarks (c-quark) that can enter in higher orders of PT through 
internal loops . 
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An important and convenient property for the analytic analysis of the above 
introduced polarization functions ^f/v/Ai^'^) absence of the so-called kine- 

matical singularities because no additional factors of momenta appear in the 
defining Eq. (|l|). For these polarization functions the dispersion relations are 
valid that describe the physical states contributing to the correlator (|l]) 

nlivM(^ ) = 1 sub (2) 

where I = 1,2, and the proper powers of s are introduced in the definition of 
R{s) to make the spectral densities positive and dimensionless. In perturbation 
theory, the threshold is at Sq = (M; + MjY (where Mj denotes the pole mass of a 
quark) while the true thresholds are, say, 4m^ for n[]] and for n|^]. It should 
be noted that the spectral density R^'^^ contains contributions from spin one as 
well as from spin zero intermediate states. The spectral density R^'^'' that is free 
from contributions of hadronic states with spin 1 is defined by 

i?(0)^^(2)_^(l)_ 

Another useful representation of the tensor Il^l^j{q) in terms of scalar functions 
reads 

n)2i (q, TUj, m, /X, a,) = (-^^.^g^ + q^q^)u''^^y/^{q^) + g;.g,.nj°^^/^(g^) (3) 

where the correlator is decomposed into the components n^°yy^(g^) that con- 
tains contributions of the states with the angular momentum J = and J = 1 
respectively. 

A direct comparison of (|l|) and @ leads us to the following relations 

n(^) = -nw /q\ n(°) = ni^i + nW/?'- (4) 

In general n[^l(0) may be different from zero which implies a kinematical singular- 
ity (pole) in both n*^^^ and 11'^'') that describe the spin structure of the correlator 



2.2 Ward identity 

The divergence of the (axial)vector current is known through equations of motion 
for the fields in QCD and is given by (pseudo) scalar two-quark operators in the 
case of massive quarks. In the massless limit both (nonsinglet) axial and vector 
currents are conserved. Nevertheless their behaviour is different and, for the axial 
current, is governed by the spontaneous violation of chiral symmetry and the 
existence of the massless excitation accompanying this violation - the Goldstone 
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boson, in our case the pion or the kaon. The (axial)vector and (pseudo) scalar 
correlators are connected through a Ward identity 

Qi^Qu'nllfjiq) = {rrii T mj^Ul^^iq) + (m^ T mj){{ipitpi) T (^jV'j)) (5) 

where 

t J rfa;e^'^^(T[ jf ^(x)(jf ^)t(O) ]) = nff (g) = Q'U^''\q) (6) 

and j^j^ = qi{i'^^)qj-, = —q^- The dimensionless function Y\!:^-^\q) is related 
to the correlator of (pseudo) scalar currents n^^^(g) by one power of Q^. 



2.3 Contour integrals and r lepton observables 

The hadronic decay rate of the r lepton is obtained by integrating the absorptive 
parts of the spectral functions with respect to the invariant hadronic mass. Cor- 
responding to two different tensor decompositions (|l|), (^ two different integral 
representations can be obtained. The first one displays the structure of hadronic 
contributions classified according to their spin 



where 



^'r dS 



1 



s 



R^^^ — \yud\^{R[idy + R[id,A) + \^us\^{R\isy + R\is,a)i ^ — 0, 1 



(7) 



(J) 



The representation of the total decay rate through the absorptive parts of the 
structure functions I!'^' and nt^^ is simpler from the point of view of its analytic 
properties for continuation into the complex plane and reads 



_ /-A^r ds / S Y 



^,R^'\s) + R^-\s) 



(9) 



Due to the analyticity of n^^'^' in the cut complex s-plane (the absence of sin- 
gularities away from the physical cut, even of kinematical singularities at the 
origin) R^ can be expressed as the contour integral along a circle C of the radius 

Isl = Ml 



Rr = Gin 



ds 

\s\=M? 



s 



(10) 



Ml 



As the behaviour of n[^'^](s) along the "the large circle" of radius |s| 
is assumed to be reliably evaluated within pQCD, representation (0) leads to 



4 



a well-defined pQCD prediction for Rr. Unfortunately, this is not true for the 
spin-separated parts. Indeed, the direct use of (H) leads to 



Qiir 



ds 



M2 



^fi-4.V 



2^ n«(^) + nW(o)/. 



:iii 



|s|=M2 M2 



M2 ; 



n(°)(s) -nW(o)/s' 



where the contribution of the singularity at the origin (proportional to n[^l(0)) 
has to be included. A nonvanishing value of n[^l(0) is certainly a nonperturbative 
constant. Thus, within pQCD we cannot predict the decay rates R[^'^^ separately. 
In the massless limit 11'^°^ = within perturbation theory and R^^^ is saturated 
by n[^l(0) corresponding to the massless pion (kaon) pole for the axial part of the 
correlator. 

On the other hand, the unknown constant drops out if one considers moments 



^(l,0)fc,«( 



■So J 



^0 ds ( \ i ' 

Mfly ~Mf; \M^, 



ds 



(12) 



with > 0, / > 1. (Note that the moments introduced in [|T5l are related to ours 

as = rW'''^ + R(o)k,i 

The decay rate Rr may be expressed as the sum of different contributions 
corresponding to Cabibbo suppressed or allowed decay modes, vector or axial 
vector contributions and the mass dimension of the corrections 



Rt — Rt,v + Rt,a + Rt,s 



(13) 



with 



R 



V 



3 



R, 



Rs 



ud 



2 

31 K J' 



1 + 5o+ ^V,ud,D 
\ D=2,4,... 

^ ' 1 + (5o + ^A,ud,D 



(14) 



_D=2,4,... 

us,D 

D=2A,... 



and 



Here D indicates the mass dimension of the fractional corrections Sy/A,ij,D, 
Sij^D denotes the average of the vector and the axial vector contributions: Sij^D = 
{^v,ij,D + SA,ij,D)/'^- If a decomposition into different spin/parity contributions 
is made or a particular pattern of moments is considered then we will use the 
corresponding obvious generalization of (0). For instance. 



R 



{l)kl 
S,V 



O-kl \ Vu. 



i + 5o'+ E 



(l),kl 
V,us,D 



(15) 
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and 

4T = \yus? ( E C'^dI • (16) 

\D=2,4,... / 

Thus, in our notation we have the relation 



3 Mass Terms in Perturbative QCD 

We would like to stress that inclusion of Cabibbo suppressed modes into the anal- 
ysis of observables related to r lepton decays gives not only an additional set of 
experimental data but open conceptually new possibilities because the massive 
piece can be measured in conjunction with the massless contribution thus provid- 
ing a strict normalization and reducing the systematic errors of the experimental 
data. In this section we therefore compute mass corrections to the moments of 
order (m^/M^) using first finite order perturbation theory and then the resumed 
perturbation theory with the evaluation of and the quark mass treated ex- 
actly (neglecting only unknown higher order corrections to the /9-function and 
the quark mass anomalous dimension). 



3.1 Finite order analysis 

Just like the perturbative predictions for the massless correlators also the quark 
mass corrections for the vector and axial correlators are identical^ for the case 
under consideration with = ^ and mj = rriu = rrid = 0. The perturbative 
prediction for the quadratic mass corrections up to order and for arbitrary 
quark masses has been presented in for the transversal piece of the correlator. 



The longitudinal piece is related to the scalar correlator through the Ward identity 



VIZ. 



l'Ki,v/A + l"Ki,v/A = l'^^v/A = ^'n^/P + m.iiss) T {uu)). (18) 

The vacuum expectation values on the r.h.s. can be understood within the frame- 
work of perturbation theory and minimal subtraction. Then formally the last term 
in (|18D is of order m^. Working only within the second order in quark masses, the 
vacuum expectation values on the right hand side of (p!8| ) can be safely discarded 
at this point. We return to them in the section for nonperturbative contributions. 
Thus, the 0{ml) contribution to the longitudinal structure function n*^*^) can be 
taken from Ref. fl^ where the massless scalar correlator has been computed at a^. 



^This is strictly true only for the perturbative contributions. 
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The resulting polarization functions Hyy^, 
in the form 



Z = 1, 2 is conveniently represented 



{Q 



usy/A 



+ 



D>2 



(19) 

Here the first term on the rhs corresponds to the massless limit while the first 
term in the sum stands for quadratic mass corrections. A similar decomposition is 
assumed for the polarization functions Il^J^\q). The results for both polarization 
functions with D = 2 read 



n 



[1] 

V,2 



+ 



2mf <i l^Q+ 



+ — 



as 

7T 



a. 



18841 



25 
T 
1 



4C(3) + 3W + ^ 



2 



3607 



-vr 



C(3) 



1265 
~27~ 



C(5) 



432 360 54 

4591. 35 22 2 17 3 



+ 



+ 

4633 
~36~ 



144 2 
1967833 1 



5184 
475 



vr 

36 



3 

11795 
24 



C(3) + ^C(5) 



108 



cm,Q + -iIq 



79,2 , 2_21 
96 



(20) 



n 



[2] 
V,2 



-Ami I 1+- 



+ 21 



+ - 



TT 





n 








n 







13981 323 
432 ^ 



520^.. 35, 
^C(5) + yW + -/ 



4 vq 



/ 14485 3659 
U [-^ + C(3) 



3380 ,,^, 1643, 
C(5) + 



27 



24 



221 2 ' 
+ 24 



(21) 



Here /^q = In^, the mass rus as well as QCD coupling constant as are 
understood to be taken at a generic value of the t' Hooft mass fi. All correlators 
are renormalized within MS-scheme. Note that terms of order are known only 



with "logarithmic" accuracy, that is the constant parts ^3^' and k^^' in the large 
Q behavior of the corresponding correlators are not available^. 

■^In fact the very calculation of the constant parts is well beyond the present calculational 
techniques. 



J2] 
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This means that we do not know the O(m^ag) contributions to Rr] however 
these constants (as well as the also unknown "low-energy" constant nW(0)) do 
not appear in the moments R^''^^^''- with A; > 0, / > 1. 

Let us discuss some concrete results (neglecting for the moment any nonper- 
turbative condensate contributions and quartic mass corrections as well). Since 
the perturbative result, in particular the terms, do not differentiate between 
vector and axial vector channels we will consider in the following their sum. The 
mass corrections to the moments of spin and spin 1 final state distributions can 
now be cast into the following form (m^ = rris^Mr), as = as{Mr)) 



In view of the large size of the coefficients bi we attempt to reduce the higher 
order corrections by adopting the value 1 GeV for the renormalization scale of 
the running mass and define an alternative set of coefficients hi through 

C^^o^^{l.S.^..(^)%.3(^)V (23, 

The results are listed in Table 1. The same expansions can also be obtained for 
the correction to the spin zero and spin one parts separately. In view of the 
additional nonperturbative contribution ~ n[^l(0) which appears in the lowest 
order moment {I = 0) due to the spin separation only the results for / > 1 
are given in Tables 2-3 for these latter cases. The additional nonperturbative 
contribution ~ n[^l(0) in the axial vector current can be estimated in the chiral 
limit and happens to be connected to the contribution of the pseudoscalar meson 
(pion or kaon) to the correlator. The distinguished role of the pion is due to 
its nature as a Goldstone particle related to spontaneous violation of the chiral 
symmetry of QCD. When the explicit violation of the chiral symmetry given by 
nonzero values of quark masses is small these masses provide small corrections 
to the massless limit (pion dominance) that can be accounted for on a regular 
basis within Chiral Perturbation Theory |T^. As for corresponding quantity in 
the vector channel there is no solid physical arguments for estimating its value 
with finite masses of quarks though it vanishes in the massless limit due to vector 
current conservation and the pattern of spontaneous symmetry breaking. 



3.2 Resummation of effects of running 

The rapid change of coefficients of perturbation theory expansions for different 
moments is caused by the running of the coupling constant and the mass along 
the contour of integration. The resummation of these effects can be performed in 
all orders of The technique for the massless case is described in the literature 
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0, |T^, |T9l so here we concentrate on the massive case. It introduces no much 
new technically but extends the freedom of choosing the resummation procedure 
due to additional parameter - the running mass of a quark. 

Let us, nevertheless, recall the central idea for the massless case. The main 
object to start with is the Adler's D-function for the transverse part of the corre- 
lator. The orders of perturbation theory in the MS scheme are formally counted 
after the renormalization group improvement of the series. The information from 
the perturbative treatment is contained in the string of the coefficients of powers 
of the running coupling. The polarization function that appears in the integral 
over the large circle is then restored from the Z)-function (including renormaliza- 
tion group improvement) by solving the corresponding evolution equation exactly 
with the j3 function taken to a given fixed order in as- Let us demonstrate the 
procedure in a simple example. Consider a D-function with its leading term nor- 
mahzed to unity (here and below we are using a = as/ii) 

D{Q^) = 1 + a(Q2) + kiaiQ^ + haiQ^f + ha{Q^f + .... (24) 

It is connected to n(Q^) by 

Dm = -^'^n(g^). (25) 



In leading order for the (3 function 



d 



/3(a) = -/3oa^ Q^^,am=P{am) (26) 

the polarization function can be completely restored in the closed form. Direct 
integration of Eq. ( ^Bj) gives the improved polarization function 

n(Q^) = in(fi7Q=) + n(,^=) 

that ought to be integrated along the contour in the complex plane with 

"^^'^ ^ i + M^h^W) ^^^^ 

which can be continued into the complex plane. The normalization point 
H is conveniently chosen to be equal to the r lepton mass /i = Mr- The part 
independent of (the integration constant 11 (/z^)) does not contribute to the 
integral. This is also obvious from the spectral integral itself because n(/i^) being 
a constant has no discontinuity across the physical cut in the complex plane. 
The first term in ( p7|) is simply the partonic contribution ln(/i^/(5^). 
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The generalization to higher orders of /3 function is straightforward. In second 
order the integration of Eq. (pSf ) can be still performed explicitly while for third 
and fourth order of the (3 function numerical integration is more convenient. 

To make things clearer let us show the explicit example in the second order 
for the (3 function 

(3{a) = -poa^ - Pia^. (29) 

(In higher orders the /3-function is scheme dependent. However, there is always 
a scheme where this is the complete expression for the (3 function - the t' Hooft 
scheme. Then the convention could be to resum always in this scheme that is 
justified by technical simplicity. Nevertheless, we however stick to the MS- scheme 
in the present paper.) A contribution to the D-function of the form AD{Q'^) = 
a{Q'^Y leads after integration of Eq. (^) to the corresponding contribution to 
the polarization function An(Q^) of the form 

An(g2) = 1 In [am + 1^ + An(/.2) (30) 

as can be easily checked by direct differentiation. The improved polarization func- 
tion An(Q^) ( |30D can then be used for integration along the contour. The remark 
about the constant terms applies again. 

The appearance of mass introduces another freedom in the choice of the basic 
quantities that accumulate the perturbative information. The actual procedure is 
described below. As basic quantities we choose nW(g^) and n[^^](g^). The renor- 
malization of the pieces proportional to ml is different for iP^ and nt^l. The 
second one n|22(Q^); is scale-invariant and the renormalization group improve- 
ment can be performed directly 

Q'n2,2(Q') = 4'^rn^(Q2)(i + ^[2]^(g2) _^ _ _ ^31) 

In contrast 11^^' is not renormalized multiplicatively and the corresponding renor- 
malization group equation is not uniform i.e it has a free term. The problem is 
solved by introducing the corresponding D-function D^^^Q"^) by one differentia- 
tion with respect to Q^. The result is 

Then we proceed as explained before in construction of the polarization operator. 
The running of the mass as taken into account through the renormalization group 
equation 

Q'^jT^msiQ'^) = 'yrnia{Q^))ms{Q'^), 7^(0) = -70a - 7ia^ -.. . (33) 
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with the solution 

rrisiQ ) =m,(/i ) exp / . (34) 

Subsequently, the integration along the contour can be performed directly. 

The explicit formula in the leading order of the /3 function for n[^l(Q^) is 
easily found. Having solved the Eq. (|3^) for the leading order /? and 7 functions 
(or using the explicit formula (@)) we have 

Q'H^uUQ') = kt^^W) { (1 + + . . .) (35) 



where a(Q^) is the solution to the renormalization group Eq. (pHf). 

For the second amplitude the result is slightly more complicated, namely 

n2.2(g') 




(36) 

These formulae should be substituted in (|10D,(|ll]) and integrated along the 
contour. The generalization to higher orders of the (3 and 7 functions is straight- 
forward. Again up to the second order expansion of these functions in ex- 
plicit analytical integration can be performed in terms of elementary functions for 
some amplitudes. (Third order also allows some integrations in terms of elemen- 
tary functions but formulas become too awkward and the numerical treatment is 
preferable.) The coefficients rff^ can be inferred from the explicit expression for 
the polarization function (^). They accumulate the whole perturbative informa- 
tion because the /9-function and the quark anomalous dimension necessary for 
the restoration of the full expressions (|20|) and ( pT]) are known. 

We perform the analysis along these lines for the same moments as in the 
finite order of perturbation theory with the available coefficients for nt^^ and for 
the D functions in the case of 11 [^1 and with (5 and 7 functions from first to fourth 



order |20|, Ell 



For the presentation it is convenient to order the results according to the 
corresponding contributions from D /Tl functions and to attach a power of ^^(M^) 
to every consequent term mark. 

For instance the typical answer analogous to Eqs. (p2D ,(p3D can then be cast 
into the form (here and below = ms{Mr)) 

~C = ^oK)^{l + + Has) (^)' + h{as) (^)'} (37) 
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where the coefficients bi{as) are obtained with resummation and become functions 
of as- These functions are not given in an exphcit form but are fixed numbers 
obtained from the actual analysis because the contour integration was not done 
in general analytic form but only numerically. The resummation procedure is 
expected to exhibit an improved behaviour of the series compared to the finite 
order analysis and to provide a definite coefficient in front of the free parameter 
ml 

The results are presented in Tables 3-6. By comparing the results of these two 
analyses one can judge the degree of improvement achieved by resummation in 
MS scheme. 



4 Non-perturbative corrections for J = mo- 
ments 

Besides the perturbative radiative corrections to Rr also nonperturbative QCD 
effects influence the hadronic r decay rate and its differential distributions ana- 
lyzed with the moments. The short distance operator product expansion (OPE) 
for the spectral functions 

Il('\s = -Q')= E 77^2 E C^'\Q',f^){0{f.)) (38) 

Z)=0,2,4, -- W ) dimO=D 

may be used to take into account both perturbative and nonperturbative contri- 
butions. Here we collect the known non-perturbative corrections to the structure 
function n^l^. 



4.1 Dimension = 4 and D = 6 Corrections 

Due to Ward identity (|18|) the = 4 part of the polarization function Byj^ can 



be cast into the following convenient form 

.2tt /r^2 



(39) 



where Il(^s/P),2 is Il25[| 



n 



S/P,2 

as 

+ — 



m' 



{ - 4 - 4/ 



+ 



33109 
108 



135 



7r^ + 



1978 



620 



C(3)-— C(5) 



12 



5065, 308 ^,^,, 



2 



50 3 
3 Vq 



-4 1. 



3.5251— + 28.0923 ( — 

TT 



(40) 
(41) 



In the last line we have set fi^ = Q"^. Note that there are no corrections to the 
coefficient function of the quark condensate term in Eq. ( ^91) because it comes 
from equal-time commutation relation for the quark fields that are just initial 
conditions for the field theory. 

It is clear from ( plSf ) that these power corrections are essentially given by the 

( S) 

corresponding D = 4 terms of the OPE for the structure function 11^^^ which are 
well-known (see, e.g. Refs. [23, EH] and references therein). As a result we have 



n 



(0) 



± 



+ 



a. 11 



vr 2 



TT 



TT 



(42) 



1 + + 2— /;,Q ) ms{uu) 

on IT 



11a 



a. 



+ -il + —^ + 2—l^Q)ms{ss) 



3 TT 



vr 



167r2 



a 



mUl- 2l^Q + 8C(3) - 6 - 4/ 



TT 



The corresponding contributions to the (fc, /) = (0, 1) moment read (to save space 
we have partially converted the coefficient in front of term to numbers and 
put /i = Mt- so that below = ^^(Mt-) and = ms(MT-)). 



5. 
and 



(0),01 
us A 



45 ml 



QJ 

1. + 4.77815— + 31.047^ 

TT 



127r' 



{{tjIsSs) =F {rrisUu)) 



rrvi 



(43) 



r(0),01 
-'us, 6 



M6 



-3 7rMl + 



a.o 11 



,a 



T 2477^ 1 + 



- 1277^ 1 + 



TT 2 

14 a 



3 TT 

11a, 



3 TT 



TT 

{nisUu) 
{rrisss) 



9 4 

-m„ 
2 ' 



1 + — 

TT 



8C(3) + %^ 



22 



(44) 



Note, please, that in Eqs. 



we have used the very quark and gluon conden- 



sates normalized at the "natural" scale fi = Mj. for our RG invariant condensates 
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(for more details, see [|, 0, [2^): 

(^GG) = (^ii^GG(^ = Mr)), 
{rrii )= {mi{Mr) ^> ^ (/i = ) ) . 

5 Discussion of numerical results 

Let us first consider the mass corrections to the total rate within fixed order per- 
turbation theory. The coefficients hi are already fairly large for all the moments. 
Even worse, however, is the fact that the series grows dramatically for all the 
moments, rendering a mass determination from the total rate or the correspond- 
ing momenta potentially unreliable. The transition from ms{Mr) to ms{l GeV) 
reduces the coefficients only marginally. From Table 2 it is evident that the large 
corrections originate mainly from the spin zero contribution to the moments, 
leaving significantly smaller QCD corrections to the spin one part. Hence, it is 
in principle possible to determine iris from the spin one part alone. Alternatively 
one might try to combine spin zero and spin one moments with different relative 
coefficients to artificially decrease the QCD corrections further. However, in the 
absence of any physically motivated reasonings we refrain from such an approach. 
The main problem of the strange quark determination and even the total pertur- 
bative analysis of the Cabibbo suppressed r lepton decays is the interpretation 
of the perturbation theory series for ml corrections. In passing we note that also 
perturbative corrections to the power suppressed terms proportional to and 



rrig (see Eqs. (|^,|^)) might provide us with another example of bad behaviour of 
higher order terms. However, since the 0{al corrections to the spin one (and, as 
a consequence, to the total R-r) are available we concentrate below at discussing 
the ml terms only. 

The problem is fairly obvious from a consideration of the two lowest order 



moments. Let us use as{Mr) = 0.334 in computing the size of the perturbative 
contributions. The order of magnitude of the unknown coefficient ^3^' in Eq. (pi]) 
is estimated on the basis of a geometric series to amount to fcf = 0±19.6V2.33 = 
± {k^2^)'^ / kf'^ = ± 164.4. For fixed order we thus find series is thus given by 
(below a = as{Mr)/ii and = ms{Mr)) 

2 

= -8^(1. + 5.33a + 46.0a2 + 284a^ + 0.75a^A;P) 

T 

2 

Tfi 

= -8—^(1. + 0.567 + 0.520 + 0.341 ± 0.148) 

2 

TT) 

= -8^(2.4 ±0.5), (46) 
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where we have assumed the (maximal!) value of the 0{a^) term as an estimate 
of the theoretical uncertainty (this convention will be used also below). For the 
"contour improved" series one obtains 

2 

= -8^(1.44 + 3.65a + 30.9a2 + 72.20=^ + 1.180^4^1) 

2 

TTl 

= -8^(1.44 + 0.389 + 0.349 + 0.0867 ± 0.234) 

2 

= -8-1(2.26 ±0.32). (47) 

T 

Comparing the "improved" with the finite order analysis one observes that higher 
orders give numerically smaller contributions although the apparent convergence 
is rather marginal as well. The ratio between two correction terms ^u^2/^us2 ~ 
0.95 shows a relatively stable behaviour. The results within resummation tech- 
nique are stable within the allowed range of as = 0.312 — 0.356 as checked by 
direct computation. Although no firm prediction can be made as a consequence 
of the large corrections, it is at least encouraging to observe that the sign of the 
mass correction remains preserved and the consecutive terms exhibit a marginal 
decrease. Eq. (^Tf ), with the uncertainty increased by perhaps a factor 2 could be 
considered as a reasonable estimate of the strange mass corrections. 

In fixed order approximation the moments with I > 1 are independent of the 
constant /C3 , a residual dependence remains, however, in the "contour improved" 
treatment 

Cs,2 = -^^(1- 4.17a -113.0^- 1820.0=^) 
= --—1(1 - 0.443- 1.27- 2.19) 



3M2 

5 ml , 

---^(-2.9 ±2.2), (4J 



and 



~C.2 = -^^(-2.26- 14.7a - 204.a2-171.a3-13.3a3A;P) 
^^(-2.26 - 1.56 - 2.3 - 0.206 T 2.62) 



3M2 
5 ml 



-6.3 ±2.8). (49) 



The (0, 1) moments thus exhibit a rapid growth of the coefficients and, at the 
same time, with S^l^2/^ul,2 = 2.2 a strong dependence on the improvement proce- 
dure. This comparison shows that there is no consistent prediction in MS scheme 
for this observable - the first moment of the differential rate. 
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Now we turn to to the contributions of the spin one and spin zero separately. 
As stated before, the lowest moments (/ = 0) of the spin-dependent functions 
depend on the nonperturbative quantity n[^l(0). For the spin one part and for 
{k, I) = (0, 1) we find 



i 

■M2' 



= -5^(1. + 4.83a + 35.7a2 + 276.a3^ 



2 

mz 



-5—^(1. + 0.514 + 0.404 + 0.331) 

Mi 



2 



-5^(2.25 ±0.33) (50) 



and 



2 

^Woi ^ _5^(i_37 + 2.55a + IG.la^ + 135a^) 

9 

TTl 

= -5^(1.37 + 0.271 + 0.182 + 0.163) 

2 

TTl 

= -5^(2.0 + 0.2). (51) 

Note that spin 1 contribution is determined by the component n^^^ alone and is 
known up to third order. Clearly, this series is decreasing in a reasonable way 
(comparable to the behaviour of 5^° 2) sjid, at the same time, only moderately 
dependent on the improvement prescription with S^]^^ /S^^]^^ = 0.89. On the basis 
of Eq. (|5T|) this moment might well serve for a reliable determination, with a 
sufficiently careful interpretation of the theoretical uncertainty. 

The corresponding spin zero part is, per se, proportional to (not counting 
the "condensate" contributions discussed in Section 4) and thus could be consid- 
ered as ideal for a measurement of m^. However, the behaviour of the perturbative 



series 



2M3 

1. + 0.992 + 1.24 + 1.59) 



2 



and 



M2 

--^(4.8 + 1.6) (52) 

2 Mr 



3 m"^ 

^uJf = 977|(3-19 + 11.2a + 126.a2 + 289.a3 + 6.63a^kf 
3 

= -—^(3.19 + 1.19 + 1.42 + 0.347+ 1.31) 
2 
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shows a rapid growth of the coefficients. The series is not expected to provide an 
accurate prediction for the mass effects. (The same rapid growth of the coefficients 
in the perturbative series is present for the scalar correlator related to r{H —>■ bb) 
[p!7[] .) Nevertheless, an interesting estimate could also be deduced from the spin 
zero contributions, in particular after resummation. For the spin zero, spin one 
and the total rate the following relation holds 

xoi _ r(i)oi 20 (0)01 

Besides the observables themselves one has also look at convergence of the 
13 and 7 functions that determine the running along the contour. In the present 
case 

P{a) = -2.250^(1 + 1.78a + 4.470^ + 21. Oa^) 

or at a = 0.1 

p{0.1) = -0.0225(1 + 0.18 + 0.045 + 0.021) 

which is quite good. For the 7 function, however, 

7(a) = -a(l + 3.79a + 12.40^ + 44.3a=^) 

and the convergence is marginally acceptable. 

The resumed series behave in general better than those of finite order. How- 
ever, for the mass corrections they still do not satisfy the heuristic criteria of 
convergence. In practice the resummation maintains the convergence pattern of 
the corresponding D-function. For the mass corrections the D-functions them- 
selves exhibit rapidly growing coefficients of the perturbative series, whence the 
resummation does not lead to a significant improvement. Numerically the con- 
vergence for D-functions in MS scheme is marginal. For the present case the 
D-function for n[^](Q^) is given by 

D^'W) = mimii + ^« + ( w - f 

0/ 1967833 vr^ 11795,,^, 33475,,^ A 

+a C(3) + C(5) 

V 5184 36 24 ^ 108 ^ 7 

= m2(Q2)(i + 1.67a + 10.840^ + 107.530^) 

whereas a much better pattern of convergence is observed for massless part 

r^-qm^U .3/58057 779 75 



DoiQ') = l + a + a'^—- 9C(3) j + a' (^^^^ ^C(3) + y C(5) 

= 1 + a + 1.64a2 + 6.37a^ . 

This is the reason why the precise determination of the strong coupling constant 
from the r lepton lifetime is possible. 
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6 Summary 



Tau decays in the AS = 1 channel are sensitive to the mass of the strange quark. 
Moments of the four spectral functions corresponding to spin zero and spin one 
contributions, induced by vector and axial vector currents can be evaluated in 
perturbative QCD. With this motivation we have evaluated the QCD corrections 
up to order for the rate and up to order for the higher moments of the 
spectral function, separated according to the spin zero and spin one contributions. 
Resummation of the RG enhanced terms has been applied to the various monets. 

At first glance, the spin zero part of the rate and of the moments appears to 
be particularly promising for an rris determination, since this piece is multiplied 
by in the parton model. However, as a consequence of the spin zero separation 
only moments with / > can be calculated in pQCD. In addition, it turns out 
that the perturbative corrections to the spin zero part are remarkably large, 
invalidating the straightforward extraction of from such an analysis. 

Also the total rate, i.e. the combination of spin zero and spin one contributions, 
exhibits a small mass dependent term which can, in principle, also be used to 
determine m^. However, again, one observes remarkably large QCD corrections 
which can be traced — at least for the higher moments — to the spin zero part. 
This is clearly visible by considering the mass terms of the separated spin one 
piece, which exhibit reasonably smaller coefficients of the first, second and the 
third order corrections. 

Resummation leads to a modest improvement of the apparent convergence 
of the series, such that the lowest moments can be predicted with acceptable 
accuracy. Assuming sufficiently precise data, a combined analysis of the spin zero 
and one moments might thus lead to a reliable determination of the strange quark 
mass. 

We do not discuss /i (renormalization scale) dependence in MS scheme as a 
special topic fixing it always at fi = M^. Our estimate of the error for "contour 
improved" analysis (2.26 ± 0.6) is rather conservative and covers all reasonable 
change of /x. The more general topic of investigating the entire scheme dependence 
of this quantity is outside the scope of the paper. 

To conclude we note that the large value for the correction factor 5°° 2/(^8) 
leads us to a reduction of the m^-value as determined from the data on 
the basis of the earlier calculation [^. The reduction by about 15% leads to 
ms(M^) = (150 ± SOexp ± 20t/i) MeV. In view of the large corrections the theo- 
retical uncertainity can just be considered as a guess based on Eq. ( ^T]) and the 
subsequent discussion. This corresponds ms(l GeV) = 200 ± 40exp ± 30(h, MeV in 
good agreement with other determinations (see, e.g. ^7\, ^ Q). 
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Table 1: corrections to Rr and its moments; for the entries with I — 
unknown coefficinet k^^'^ is set to zero 
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Table 2: corrections to the spin zero part of Rr 
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Table 3: corrections to the spin one part of Rr 
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Table 4: corrections to Rr and its moments; q;s(Mt-) = 0.334; the unknown 
coefficient kf^ is set to zero. 
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Table 5: corrections to R^^ and its moments; q;s(Mt-) = 0.334, the unknown 
coefficient ^3^^ is set to zero. 
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Table 6: corrections to R^^ and its moments; Q;s(Mr) = 0.334. 
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